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Abstract: We present an efficient method for computing the zero frequency limit of trans- 
port coefficients in strongly coupled field theories described holographically by higher deriva- 
tive gravity theories. Hydrodynamic parameters such as shear viscosity and conductivity can 
be obtained by computing residues of poles of the off-shell lagrangian density. We clarify 
in which sense these coefficients can be thought of as effective couplings at the horizon, and 
present analytic, Wald-like formulae for the shear viscosity and conductivity in a large class 
of general higher derivative lagrangians. We show how to apply our methods to systems at 
zero temperature but finite chemical potential. Our results imply that such theories satisfy 
r]/s = I/Att universally in the Einstein-Maxwell sector. Likewise, the zero frequency limit of 
the real part of the conductivity for such systems is shown to be universally zero, and we 
conjecture that higher derivative corrections in this sector do not modify this result to all 
orders in perturbation theory. 
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1. Introduction 



In recent years, with the AdS/CFT correspondence QlQ] extensively tested, there has been 
increasing interest in applying it as a tool to study strongly coupled systems. The case where 
the dual field theory has conformal symmetry is of particular interest as it could provide 
information on real-world systems such as the quark-gluon plasma and materials at quantum 
critical points. For recent reviews and further references on these topics see for instance 

iii- 

According to the AdS/CFT dictionary, a conformal field theory at finite temperature is 
dual to a black hole living in an asymptotically AdS space. The thermodynamics of the dual 
field theory then correspond to the thermodynamics of the black hole. It can also be shown 
that small hydrodynamic fluctuations around thermodynamical equilibrium map onto small 
perturbations of the gravitational background in a precise fashion Q . 

A hydrodynamic description is relevant at wavelengths and frequencies that are large 
compared to some typical microscopic scale of the theory. In this approach one describes the 
system by a set of transport coefficients. Determining these quantities from first principles 
is in general quite difficult. Typically we are interested in transport of conserved charges 
in the field theory, and in this case Kubo formulas can be used to read off the coefficients 
from low frequency poles in the thermal retarded Green's functions of conserved currents. 
This requires a study of real-time physics at finite temperature. While this can be achieved 
with some difficulty at weak coupling, at strong coupling the situation is thornier, as lattice 
methods are intrinsically euclidean. For some recent results on the shear viscosity of Yang- 
Mills theory plasmas see for instance ||^]. 

In contrast, in a holographic context these computations turn out to be remarkably 
simple. The prescription for calculating thermal retarded Green's functions was first set out by 
[IIIQ, and later more formally justified in |[TT[|,|p^. Using this approach, all hydrodynamic 
transport coefficients up to second order have been computed in a large class of theories, as 
well as their first string theoretical corrections @|g[|l3|[|l§0|l8|[|l9||O|. 



In a beautiful paper, the authors of |21] showed that there is a very close link between the 
transport properties of the black hole "stretched horizon" in the membrane paradigm p^ , 
and those of the dual field theory. In this picture, the Green's function is essentially given 
by the canonical momentum associated to the bulk perturbation that sources the boundary 
conserved current. While the value of the canonical momentum at the horizon is fixed by 
regularity, generically there is a radial flow to the boundary. In the hydrodynamic limit this 
flow can be trivial, and in this case the properties of the field theory plasma coincide with 
those of the membrane. This is exactly what happens for the shear viscosity, from which the 
celebrated result rj/ s = l/An can be understood as emerging from the universal properties 



of black hole horizons. Universal properties of correlators in the hydrodynamic limit had been 



noticed already in \24\. 



Extension of these methods for the computation of the shear viscosity in higher derivative 



theories have been given in |25],|26]. The approaches of those papers are somewhat comple- 
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mentary. Whereas the first focuses on the generalization of the canonical momentum method, 
the second proposes to rewrite the effective action as a two derivative theory. The problem 
with the latter is that the overall coefficient of the effective action has to be fixed by hand, 
and it is not clear what its value should be in a general higher derivative theory (e.g. with 
derivatives of curvatures). Even though the two approaches should be in principle equiva- 
lent, both procedures rely on massaging the effective action into a desired form, which for 
general higher derivative theories can become quite cumbersome. Furthermore the treatment 
of possible contributions from boundary terms is not done in a systematic fashion. 

In both these methods it is clear that generically transport coefficients associated to 
massless modes are given by effective couplings at the horizon. In this paper we show that 
there is a quick, efficient procedure for extracting these couplings from the lagrangian density. 
By evaluating the lagrangian on an off-shell perturbation, it will generically develop a pole at 
the horizon. The residue of the pole is precisely the desired transport coefficient up to a known 
factor. Boundary terms cannot contribute as in general they yield higher order singularities. 
The single pole behaviour can be traced back to horizon regularity of the generalized canonical 
momentum in the hydrodynamical approximation. This procedure for computing transport 
coefhcients, which we name the pole method, works for any higher derivative theory and 
reproduces all the previously known results in the literature in a simple fashion. 

The idea of effective horizon couplings giving transport properties was pioneered in 
[p^] , [p^j , where a specific covariant formula for the horizon graviton coupling was proposed, 
and therefore for the shear viscosity. However, it has been found in |26| that this formula 
does not yield the correct result for the well known correction |13]. In this paper we 
clarify this issue, by carefully finding an analytic expression for membrane couplings, find- 
ing disagreement with the proposal of [27|. Specifically we find formulae for both the shear 
viscosity and DC conductivity in an uncharged background for a wide class of lagrangians. 
These expressions match all previous calculations in the literature. 

Recently there has been interest in studying the transport properties for backgrounds 



at zero temperature but finite chemical potential [29|,[30|,[pl|. The relevant backgrounds 
describe a flow from the boundary field theory to an IR CFT, itself described by an AdS2 
factor in the near-horizon geometry . We show that the shear viscosity to entropy density 
ratio may be computed solely within the near horizon geometry, and therefore in Einstein 
gravity we deduce the universality of = l/47r at extremality. We show how to apply the 
pole method in this well as giving a simplified analytic formula for the shear viscosity. 

The computation of the conductivity is more interesting, as it exhibits a flow brought 
about by the background charge. As shown in |3^], in the IR CFT the real part of the 
conductivity scales like a;^ in the hydrodynamic limit. We give evidence that this result 
is actually valid generically for any higher derivative theory. The RG flow from the IR 
to the boundary theory cannot affect this result, which leads us to conclude that the low 
frequency limit of the conductivity of the dual theory is zero to all orders in the higher 
derivative expansion in the gravity-gauge sector. This is equivalently described as the non- 
renormalization of the conformal dimension of the specific scalar in IR CFT associated to 
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charge transport. 

The structure of this paper is as follows. In section 2 we describe the basic setup of 
our calculations for field theories at finite temperature. The following section clarifies the 
notion of effective horizon couplings and shows that the generalized canonical momentum 
method is generally applicable. Section 4 introduces our pole method for computing transport 
coefficients and some practical examples. In section 5 we derive analytic formulae for the 
shear viscosity and conductivity for a vast class of theories and compare them with previous 
proposals in the literature. Section 6 is devoted to the study of systems at zero temperature 
but finite chemical potential. We show that our methods are still applicable there and as 
an immediate corollary prove universality of r]/s = 1/4tt in two derivative theories. We also 
study the low frequency limit of the real part of the conductivity in general higher derivative 
theories and find evidence that it is zero to all orders in perturbation theory. We finish this 
paper with a discussion and future prospects. Some technical points are presented in the 
appendices. 

2. General Setup 

We are interested in computing transport coefficients for field theories at finite temperature 
living in d flat dimensions. We assume that, in what concerns thermal properties and hydro- 
dynamics, these theories admit an effective holographic description at strong coupling by a 
dual {d + l)-dimensional gravitational background of the form: 

ds'^ = gabdx'^dx'' = — e^^^^^dz^ + g^^dx^dx" 

g^^ = -ze'^f^'Ut'^ + e'^P'^'Ux'dxi. (2.1) 

The functions /, g, p are regular everywhere in the bulk, which means that there is a horizon 
located at z = 0. Henceforth we will denote quantities evaluated at the horizon by an index 
'0'. The field theory lives at the boundary of the space aX, z =1 and is at temperature^ 

r=^exp(/o-5o). (2.2) 

It will be convenient to define the volume of the horizon as given hy V = e^'^~^^^° . 

We shall consider perturbations of the above background described by a massless scalar 
field (l){z,x'^). According to gauge-gravity duality, such a field sources an operator O in the 
dual field theory via a coupling / d'^xO{x^)(l)B{x^)^ where (f)B denotes the boundary value of 
the perturbation. 

Prom the field theory perspective, we have simply some source coupling to the operator 
O. For small enough fields, and going to momentum space, linearized theory then determines 
the expectation value of O as 

{0{k)) = GR{k)Mk)- (2.3) 

^Tho temperature can bo obtained by fixing the periodicity of the Wick-rotated time coordinate as to 
eliminate the conical singularity at a = 0. 
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In the low frequency limit and at zero spatial momentum one expects {0{uj)) = i^^uicj). There- 
fore, to obtain ^ we can use the Kubo formula 

^ = lim ilmGij(Lj,k = 0). (2.4) 

Finding the transport coefficient ^ then amounts to computing the retarded Green's function 
in the small frequency limit and at zero spatial momentum. At weak coupling it is possible 
to directly compute this Green's function in the field theory. However, at strong coupling it 
is simpler to perform this computation holographically. 



3. Retarded Green's functions in the hydrodynamic approximation 

In this section we show how to holographically compute thermal retarded Green's functions. 
We start by briefly reviewing the canonical momentum method of and then show how to 
generalize it for higher derivative theories, including a careful treatment of possible bound- 
ary terms. We study the regularity conditions at the horizon, and show precisely how the 
transport coefficient ^ is related to an effective coupling at the horizon, along the lines of |27]. 

3.1 Two derivative theory 

The procedure that must be taken to compute retarded Green's functions via gauge-gravity 
duality is by now relatively well understood. Consider adding to the background ( p.l[ ) a small 
perturbation described by an action: 



ST = -2 J d'xdz^{V<t>r. (3.1) 

We now go to momentum space, and set the spatial momentum to zero. This is accomplished 
by choosing the ansatz 

^ ' ' ''^"^^ (3.2) 



(f>{z,x) = J ■^-:^(f>^{z)6{k)e' 
with = (— a',k). The action becomes 



Sf = /n^^V^ [\z^f^{g^^^U^)^'-.(.^)+g'WM^)^-M)) , (3.3) 



Following [p^], the retarded Green's function is given by the formula 



Gnik) = hm ^ (3.4) 



where 'Kuj{z) is the radial canonical momentum associated with (t)w{z): 



5S 



(2) 
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Now consider the small frequency limit where we take w/T — > keeping TTujiz) fixed. In these 
conditions the radial evolution of both canonical momentum and uxh is trivial lEltl (as we will 



show in the next section), and so we are free to evaluate the ratio in (3.4) at any radius we 



wish. We will see that it is convenient to choose the horizon for this so that, 

e=lim^, (3.6) 

where we have used uj(j){z) = iV(f)Q. 

3.2 The generalized canonical momentum 

We would like to generalize this procedure to theories containing higher derivatives. Such 
theories can arise for instance by considering string theory corrections to the supergravity 
action. It is implicit that there should be some small parameters controlling the size of such 
corrections. Then the higher derivative terms are perturbatively small, and the equations 
of motion can always be recast as two derivative equations. In this paper we shall always 
assume that whenever there higher derivative terms present, they are always parameterically 
smaller than the two derivative terms. 

We will continue to restrict ourselves to massless fields, and work in Fourier space at zero 
spatial momenta. In this case, the most general quadratic action can always be written as: 

i=l 

The first piece in the above is the "radial" action. This contains all the terms in the action 
without any time derivatives, or equivalently the action that one obtains by setting a; = 0. 
The second piece contains the remaining terms, and by time reversal invariance it is necessarily 
proportional to o;^. Finally Sb contains boundary terms necessary to make the variational 
problem well defined^. We will consider radial actions of the form 



Si.)= I dz\ ^„,^(z),^(r+^)(.)0Lr^^(z) 1 , (3.^ 

where (j)u\z) indicates {dzY4'uj{z). We should emphasize that in picking the specific form 
of the action (^), (|3.8D , the only constraint is that the perturbation cj) should be massless. 
Under these conditions, it is always possible to integrate by parts so that we get an action of 
this form. We will see that in applications these manipulations will be unneccessary. 




For an action given by (3.7),( |3.8| ) the boundary action necessarily contains only three 
types of terms. The first is simply Bq{z)(IP'. This type of boundary terms does not contribute 
to the imaginary part of the Green's function, so we will not be concerned with them. Then 



^For higher derivative theories, this has to be done perturbatively in the parameters controlling the higher 



derivative terms. See e.g. jl^ for an example on how this is done in practice. 
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there are boundary terms proportional to cj^, which in the small frequency limit do not 
contribute. Finally there are terms of the form i3„ These terms must 

always have at least one derivative on each (j)uj{z). To see this, consider for instance the 
following term in the radial action, 

A,,2{zmz)^'LJz) (3.9) 

Variation leads to boundary terms 

{Ai,2{z)6(t>UzWLjz) + Ai,2{z)(t>Uz)H'-Jz)) (3.10) 

We may set the first term to zero, but we must have a boundary term of the type (jy'cj)' to 
cancel the second term. We see that no undifferentiated <j) appears. It is clear that considering 
terms in the radial action with more derivatives will lead to similar results. To sum up our 
description of the boundary action, we conclude that, in the low frequency limit, the only 
possible relevant boundary terms are such that they only contain differentiated fields. 
Under these conditions we define the generalized canonical momentum: 

n.(z) ^ — ^ (3.11) 

To apply this definition, one can simply treat dz(j)~.uj as a new field ip{z) and then treat (p, tp as 
independent fields. The canonical momentum is then the functional derivative with respect to 
the ipiz) field. There is no ambiguity in this definition, as (^(^{z) always appears differentiated 
in the radial action. With this definition in hand, we can integrate by parts so that the radial 
action becomes ^ 

Si.) = j^dz Qn^(z)0'_,(^)) (3.12) 

This rewriting will produce some extra boundary terms which are covered by our previous 
discussion. Then it is clear that the equation of motion derived from the total action (|3.7| ) is 
of the form: 

5,n^(z) =c^2^(z, </),(/-',...). (3.13) 
The Green's function is once again given by the value of the on-shell action: 

Gr(u;) = lim ^\ \ + Boundary terms. (3-14) 

z~*l (j)^{z) 

In the small frequency limit, the equation of motion implies dz^uj — 0. On the other hand, 
notice that at strictly a; = 0, a constant field (^{z) = (j)^ must be a solution, and therefore 
dz<P{z) = 0{lo). In this way we have shown that just as in the two derivative case we have 



a trivial flow from horizon to boundary, and so we are free to evaluate the ratio in (3.14) at 
the horizon as before. 

One needs to consider as well the possible contributions of boundary terms. As discussed 
above, a typical term is of the form dz[B(z)(j)'^~^^(z)(j)'^'^^(z)]. This form can be traced back to 
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the way in which we picked our radial action. But any differentiated (f){z) is necessarily 0{uj) 
as we've seen. Then automatically these boundary terms are all of order w^, and therefore 
cannot contribute to the imaginary part of the Green's function in the low frequency limit. 
Using (U) and (|3ll ) we conclude that 



«=limHliM. (3.15) 

This shows that for a large class of higher derivative theories, the imaginary part of the 
retarded Green's functions at zero spatial momentum and small frequency is still given by the 
canonical momentum term. This generalizes the results of to arbitrary higher derivative 
theories and makes clear that there can be no contributions from boundary terms. 

3.3 Near horizon behaviour 

Our results show that the transport coefficient ^ is completely determined by the horizon 
behaviour of the bulk field z). Following |^l|], we notice that infalling observers must see 
a regular field 4'{t, z) at the horizon. This means that for these observers, (j) can only depend 
on the ingoing Eddington-Finkelstein coordinate v defined by 



dv = dt+ J-—dz. (3.16) 
V 9tt 

Actually time reversal symmetry implies that another solution is possible. Namely, it could 
also be the case that 4>{t, z) depends on the outgoing Eddington-Finkelstein coordinate u. 



du = dt- J-—dz. (3.17) 

9tt 



This means that 



d.cPo = ±J-'—dt^o = Ti^^. (3.18) 
gtt 47rT z 



By composing the two possible behaviours ( p. 18 ), it must be that the equation of motion at 
the horizon is given by 

z [Airl 



Indeed this is what follows from the action ( |3.3| ) once we take the z ^ limit. This argument 
is completely independent of the specific action one is considering, and so any action for the 
perturbation should lead to the same horizon behaviour. Regularity at the horizon is the only 
constraint. If we further demand on physical grounds that the perturbation should actually 
be infalling, the unique solution to the above equation is given by 

4>uj{z) = (f>o exp (-i-^ log z) . (3.20) 



AirT 

This completely fixes the near horizon behaviour of the perturbation. 



- 8 - 



We see then that quite generally, close to the horizon the scalar perturbation must satisfy 
the universal equation (3.19) regardless of the theory one is considering. This seems to be 
related to the conformal symmetry associated to the black hole horizon. If one looks at the 
form of the action discussed in subsection 3.2 one would have no reason to believe that the 
implied equation of motion would reduce to ( 3.1S| ) in the near-horizon limit. However, one 
must remember that the various functions involved are not arbitrary, and are only effectively 
deduced from an original scalar lagrangean density, which is necessarily smooth at the horizon. 

Regularity then imposes that close to z = the effective action must be of the form 

= /n^^V^ [\z^^ {g'^(t>Uz)ct>LM) + 9''^'Mz)<l^-Uz)) , (3.21) 
J i=i J Jo ^ 

plus boundary terms. It is important to note that this form of the action does not mean 
that higher derivative terms simply drop out at the horizon. Indeed these play a crucial role 
in determining the effective horizon coupling constant k. For obtaining this two derivative 
action from the general higher derivative one could follow the approach of ||2^. For instance, 
for a single higher derivative term of coefficient 7, one applies the lowest order equation of 
motion on the 7 terms, reducing them to terms with single derivatives or less. In this way 
one obtains a second order equation which can be deduced from an action of the form above 
in the near horizon limit. The problem with this procedure is that the overall coefficient of 
the action has to be fixed by hand by comparison to the original action^. 

We will see that in practice one doesn't need to work with the equations of motion at 
all. It suffices to know that at the horizon the action must necessarily take the form ( p. 21 ). 
For such an action, the canonical momentum associated is easily computed. At the horizon, 
using ( 3.20 ) we obtain^: 

U^{z) = iuj^ — . (3.22) 

K Z 

Using ( p.l5| ) this gives 

^ = ^. (3.23) 

K 

In this way we have shown precisely the direct link between the transport coefficient ^ and 
the effective coupling at the horizon k, and how this is valid even for general higher derivative 
theories. 



4. Membrane coupling as the pole at the horizon 

By now it is clear that we are only interested in the behaviour of quantities at the horizon. 
In this section we show that there is an efficient way of directly obtaining the value of the 
canonical momentum at the horizon, and therefore the membrane or horizon coupling k. This 
is done by evaluating the lagrangian off-shell in a specific fashion. We then show how to apply 
this method to compute the shear viscosity and conductivity in higher derivative theories. 

^The authors of do this for higher derivative theories containing curvatures, but not their derivatives. 
''Note that at the horizon we have yf^g"'' T/(47rT). 
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4.1 The pole method 

Let us start by considering the near horizon form of the action, equation ( ^.21 ). Consider 



plugging into the action a perturbation of the form 

(pcuiz) = (poewi-ialogz) . (4.1) 
In the near horizon limit we obtain 

i=l ^ ' 

It is clear that on-shell the lagrangian is regular at the horizon, and indeed it is zero. This 
is because it reduces to the boundary term 9^(n^(z)(/)_^(z)), and the quantity being differ- 
entiated is a constant in the near horizon limit. However, we see that as long as we keep the 
perturbation off-shell, the lagrangian has a simple pole, the residue of which corresponds to 
the desired membrane coupling k. In particular, starting from the full action, we might have 
considered a perturbation that only depended on effectively setting to zero the LtP' term 
in the above; or alternatively, a purely time dependent perturbation, which would set the o? 
term above to zero. Either way, we can find out k from the residue at the pole, and use that 
information to get ^ via ( 3.23| ). To be precise, let us start off with an effective action for (\) of 



the form: 

Sf = J d''xdzCl'>{d,(t),dt(t)). (4.3) 

Once again our only requirement is that the lagrangian above does not depend on (p itself, 
but only on its derivatives, of which there can be an arbitrary number. Then we reach the 
following simple formulae for the transport coefficient ^: 

f = 8ttT lim Radial formula (4.4) 

Res2=o J--^? — .-...t 



$ = — SvrT lim Time formula. (4.5) 



We call this the pole method for deriving ^. There is in fact an infinity family of formulae all 
exploiting the simple pole, but the above are the two simplest. 

This method works because essentially, the canonical momentum term in the lagrangian 
is the only one which develops a simple pole. In particular, all boundary terms have higher 
order divergences. This follows from the discussion of boundary terms in subsection (|3.1| ). 
Since the action only contains differentiated fields, so do the boundary terms. Necessarily 
these will then be 0{l/z'^) at the horizon. 

On the other hand, the canonical momentum term necessarily has a simple pole, since 
'^Luiz) = 11(^(0) in the low frequency limit. Using ( |3.1^ ) we obtain 

limn^(z)(/._^(z) = -— . (4.6) 

z^o 47rJ z 
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We see that the residue of the simple pole yields the value of the canonical momentum, which 
ultimately leads to ^. Let us stress that the formulae (|4.4| ),( |45| ) are generic, and apply to any 
higher derivative lagrangian. In applications, this means for instance that we may compute 
the shear viscosity in theories with an arbitrary number of derivatives of curvatures. 

The only ingredients that go into these formulae are horizon regularity, and the require- 
ment of a massless field. While the first is generically satisfied as long as we start from a 
covariant lagrangian, the second requirement is more restrictive. Notice also that knowledge 
of boundary terms is completely irrelevant, as long as the action only contains derivatives of 
the field. This can always be achieved by doing some by parts integration on the original 
action. In practice it is worth choosing our perturbations in a way that they can only appear 
differentiated, so as to avoid this extra work. We will see this is easily done for the shear 
viscosity and conductivity in the next section. 

4.2 Application to shear viscosity. 

The shear viscosity is computed on the field theory side by computing a specific two point 
function of the energy momentum tensor. More concretely, defining the retarded propagator 

G'^'^'^uj) = -i [ dte{t){T^y{t)T^y{o))e~'"^\ (4.7) 



the shear viscosity is given by the Kubo formula 

, ImG"(c^) , , 

T7 = lim ^ —. 4.8) 

To compute the Green's function holographically we need to turn on a graviton perturbation 
hxy We focus on the shear mode channel at zero momentum, corresponding to a change of 
basis: 

dx2 dx2 + Am{x'^)dx"', (4.9) 

where m runs through every coordinate except X2 ■ In radial gauge and at zero momentum, the 
{t, z) = (f>(t, z) component decouples from all others, so here we simply set them to zero. It 

is clear from gauge invariance that (j){z) can never appear undifferentiated in the lagrangian, 

and so our formulae can be used. In particular, we need not worry about boundary terms. 
Let us see how to use the pole method to compute the shear viscosity. Start by considering 

the following action 

1 /• , /„ 12 



I d^'x^ + ^ + iL^^aRbcde^^R''''''^ (4.10) 



levrGiv 

where 7^1. The AdS-Schwarzschild metric. 
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with i = 1, .., 3 extremizes the action when 7 = 0. The associated temperature is easily found 
to be T = ro/{L^ir). In the presence of the higher derivative term this is no longer a solution 
and receives an order 7 correction. 

However, one does not need to compute the order 7 correction to the background metric 
if one is only interested in the r]/s ratio. This is a point which although known by experts, 
has not been made clear in the literature. If one is working to linear order in 7, then 
clearly it is sufficient to evaluate the higher derivative terms to the lagrangian on the lowest 



order background. On the other hand, in Einstein theory universality |23| guarantees that 
r]/s = I/Att regardless of the details of the background we're considering. 

Perturbing the metric by Ax-^ (t, z) = (l){z) , and evaluating the lagragian to quadratic 
order gives: 

^i'^ = "32^ + + 

+ i^^L'V-. + i^#<A". + F<p(J^UZ) (4.12) 



where the functions A, B, F are given in appendix B. Let us start by computing the shear 
viscosity using the canonical momentum method. Applying definition ( 3.11| ) gives 



nUz) = Act>Uz) - {BUz))' + {E<t^'L{z)Y (4.13) 

with 

A = A--B' + -D" 
2 2 

B = C --E' - D. 
2 



Plugging in the near horizon solution (3.20) gives 



n.(0) = i^ ^g^^°3^^ (l - 10247)</)o + 0(a;2), (4.14) 



from which we can read off the shear viscosity, 

Now, let us try the pole method. We simply evaluate the lagrangian density associated 
to the action ( [4.10| ) on the perturbed background metric to quadratic order in the field cj). 
Substituting (j){z) = ^^'^^^/(^''o) and taking 2; ^ we obtain 

Had we evaluated the lagrangian on a metric perturbed by Ax^{t, z) = (j){t) = e"*"^* we would 
get the same result but with an opposite sign. From this expression we read off exactly the 
same value for r] as by the canonical momentum method. Computationally this is highly 
efficient: one simply evaluates the full higher derivative lagrangian on a perturbed metric and 
looks for the pole at z = 0. 
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4.3 Application to DC conductivity 



The computation of the DC conductivity can also be achieved using our method. Starting 
from the retarded propagator, 



-iujt 



the conductivity is given by the Kubo formula 



ui^O lUJ 



(4.17) 



(4.18) 



Here e is the coupling of the weakly gauged boundary field theory [32| 
Consider a general action of the form 



S 



(4.19) 



We can equivalently define the tensor X by 



X 



abed 



g 6Fab6Fcd 



(4.20) 



We assume there is no background gauge field present. Turning on a small perturbation 
Axj^{t,z) = (l){t,z), the resulting quadratic effective action does not involve undifferentiated 
perturbations: 



S 



(2) 



25^ 



^ / d'^xdz^g-- {X'%^^y^d,^d,cP + X'-\^^g''dtcpdt^) (4.21) 
d+l J 



Applying the formulae ( ^.7D ,( |5^ it is straightforward to find: 



9: 



d—3 -y-zxi 

XX ZXl 



z=0 



9x 



'A' 



tx-\ 



(4.22) 



2 = 



where the last equality is required by regularity of the effective action at the horizon. Now 
consider the special case 

Xabcd = Qacgbd — QbcQad — §7 Cabcd (4.23) 



where Cabcd is the Weyl tensor. For the AdS-Schwarzschild metric in d + 1-dimensions we 
obtain 

{^-z? _J_ 
L2 L2 

and therefore 



^ tXl ^ ZXl 



(4.24) 
(4.25) 



in agreement with 
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5. Wald-like formulae for transport coefficients 



From our results in the previous sections, it is clear that there exists a class of transport 
coefficients which are holographically given by effective couplings at the horizon, or membrane 
couplings. There is a specific proposal in the literature for a formula that gives this coupling 



for the special case of the shear viscosity [27|. It has been shown that this formula is incomplete 
as it does not yield the correct answer for well known case of a correction [|2^ . It is therefore 
of interest to find such a formula, which we do so in this section for a large class of lagrangians. 



We also give an analogous expression for DC conductivity. Equations (5^), (|5.8| ) and ( 5.11 ) 
are the main results of this section. 

5.1 Deduction of formulae 



We start in the background {2A) and add a shear mode perturbation 

dx2 dx2 + A^{x'^)dx'^, (5.1) 



Rmnpq — R-mnpq . 6 ^ PiF^mn-l^pq] 



as before. The indices m,n,p, ... do not include X2. Under such an addition the curvature 
transforms as |34] 

3 

V 

R'mnpy — 2^ ''-^mm)] (^•^) 

where hatted quantities refer to the unperturbed background and where P acting on a tensor 
gives it the relevant symmetries for the Riemann tensor. For instance 

P[Vp(e3^F„„)] = ^ [2Vp{e^PFmn) + V„(e3^F„p) - V^(e3''F„p)] (5.3) 

We consider a general lagrangian which does not involve derivatives of curvatures. We assume 
an action of the form 

5 = -Y^^ / d^^^9C{Rabcd,F^£M'\.:) (5.4) 

The action may depend on an arbitrary number of other fields such as abelian gauge fieds and 
scalars. All we require is that the shear mode appears only through the curvature, which in 
practice requires that the field strengths do not appear differentiated. Scalar fields can have 
single derivatives on them, or an arbitrary number of box operators. 

Using the expressions above, the quadratic effective action S^"^^ for the shear mode is 
given by 



_|_2 -^amnpjbqrs 



(5.5) 
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where 



X 



abed 



Yabcd,efgh 



6X 



abed 



(5.6) 



V—g SRabed ^^efgh 

and = ^ ^ other components zero. Notice that Y^X inherit all the symmetries of the 
Riemann tensor. 

Taking Ax2{t,z) = (j){z) = and expanding the effective lagrangian near the 

horizon, we can read off the residue at the pole and obtain the shear viscosity. Our calculations 
are given in detail in appendix A. The final result is 



V 



167rGiv 



zy xy 



4 



(5.7) 



z=0 



where 



.zt 



a 



yxz ,yz 



zy zx 



Y 



xz ,yt 
zy tx 



and we have renamed xi = x,X2 = y. 



Alternatively one may take Ay{t,z) = <j){t) 



-iujt 



and reading off the residue at the 



pole get the formula, 
V 



V 



where 



2X\ - x^y,^ 



4e-29o 



(5.8) 



2=0 



~zt _ yxt ,yz yxt ,y z 
^ zy tx ' zy t x' 



These formulae are similar to Wald's formula for the entropy [^]. In particular, in the 
background and in our notation the entropy density is given by 



V 



^Gn 



X 



zt 



zt 



(5.9) 



z=0 



Notice however that in this case the double Riemann derivative type coefficients do not 
contribute. This is because the viscosity is related to metric fluctuations which involve dif- 
ferentiating the action twice with respect to the metric. 

We can now do the same sort of calculation for the conductivity. We consider an action 
for the gauge field of the form 



Sa = J d'^+'xV^ (^Ax-bed p^^p^^ ^ ^ 



(5.10) 



and once again assume there is no background gauge field turned on. Then it is straightfor- 
ward to generalize the result ( [4.22 ): 



a 



r,d—3 I 



2 ilxx 
Sd+1 



2 



dzP"' - (e^^^L^/? +{d- 1)9, 



(5.11) 



2=0 
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with 



■zt 



zy zy 



yZ ,t 

zy ty 



Just like for the shear viscosity there is also another formula obtained by taking a time 
dependent perturbation. We have tested these formulae against the generalized canonical 
momentum and pole methods and have found agreement. 

5.2 Comparison with previous proposal 

The formulae deduced in the previous sections are perturbatively valid in the coefficients 
controlling the size of higher derivative corrections, but they are correct to any desired order 
one wishes to consider. In theories containing only two derivatives, such as Lovelock gravity, 
the formulae are non-perturbatively correct. 



Our formulae for the shear viscosity are distinctively different from the proposal of |27]. 
For backgrounds of the form and in our notation the proposal is 



7/ 



V 

v-xy 



(5.12) 



z=0 



The authors of |27] found that this formula gives the correct shear viscosity for Gauss-Bonnet 
theory. Let us see precisely why this happens. Consider then Gauss-Bonnet theory in five 
dimensions. The action is given by 



with correction W 
solution 



S 



(R 



1 



IGttGn 



(5.13) 



abed ~ '^^ab + • "^^^ action is extremized by the black hole 

du^ 



-| {-N' f{u)dt' + dx'dxi) + 



with 



1 



1 + Vl - 4A 



This background can be put in the general form ( |2.1| ) by doing the substitution u 
Computing the various quantities in the formula (|5.7D we obtain 



(5.14) 

(5.15) 
1 - z. 



xy 



1-4A-32A^ 



j^zy 



zy 



zt 



z=0 



z=0 



1 - 8A 



0, 



zt 



a 



XL' 



z=0 



R 



2 = 



2 = 



20 

l2 



(5.16) 



It is immediately apparent that the formula of |27] cannot be correct, as it leads to quadratic 
terms in A. We also see that by an apparent accident, the formula gives the right result 



perturbatively in A. Plugging in these values into (5/7), quadratic terms cancel and we get 

1 /rnN 3 



(5.17) 
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which indeed is the correct value in five dimensions ||36|]||3^. To obtain this result it was 
crucial to include the double-Riemann derivative type coefficients. 

Now consider applying the formulae for the well known corrections. In this case 



(5.18) 



The computation of the various coefficients involved is done straightforwardly with the aid of 
a computer^. We obtain 



zy 



1 - 207, x'^y^ 



xy 



dza 



zt 



z=0 



z=0 



I + 2O7 



z=0 



zt 



z=0 



2 ' 



(5.19) 



Naively applying formula ( |5.12| ) gives 



1 



levrGAT \L 



(I + 2O7). 



This result had already been reported in ||2^. Our formula however gives 



(5.20) 



(5.21) 



which matches the previous calculation in the literature [13|[37|. 

This disagreement can be ultimately traced to an incorrect formula for the kinetic term for 
gravitons. It seems that this is not limited to the h% gravitons, and so in the same fashion one 
expects that it is not correct to think of the Noether charge entropy as the effective coupling 



at the horizon of h^t gravitons |28]. Our results show that even in two-derivative theories the 
graviton kinetic term receives contributions double-Riemann derivative type terms. 



6. Extremal backgrounds 

In this section we will be interested in studying the holographic duals of field theories at zero 
temperature but finite chemical potential. In this regime, the holographic dual gravitational 
background is expected to contain an extremal black hole. It is then possible to take a scaling 
limit where one focuses on the near horizon geometry |2^. We will consider theories where 
this near horizon geometry is given by AdS2 x R'^~^, supported by the flux of an abelian 
gauge field Fab- In terms of the dual field theory, the full geometry describes an RG fiow 
between the boundary d-dimensional field theory and a conformal fixed point located in the 

^Notebooks are available from the author upon request. 
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IR, characterized by the AdS2 factor. The near horizon geometry is described by the fohowing 
metric and gauge field: 

ds'^ = -VI ^-z^dr^ + + V2{dx'^). 

F,r = Q. (6.1) 

The equations of motion fix Q in terms of the AdS2 radius vi. In this background, a massless 
scalar field (l){z, r) = (/)(z)e*^'^ satisfies the equation 

riz) + ^<P'iz) + "4^^) = 0- (6.2) 

This is exactly the equation of motion one would expect close to a double pole horizon. Indeed 
for a metric of the form, 

ds^ = ^ e'3iz)dz^ + (^_,2 e2/(.)^^2 ^ e^Piz)dxidx^') , (6.3) 

regularity in Eddington-Finkelstein coordinates at the horizon imposes 



Mo = ±,r^d,cPo = T-% (6.4) 
V 9tt 1^ 

where we have defined the chemical potential = e-^'^^^'^/L. By composing the two possible 
behaviours we get back ( |6.2| ) after a suitable redefinition of oj. The solution to ( |6.2| ) is given 
by 

</.(z) = (/>oexp(^±^y (6.5) 

From the point of view of the full geometry, this only fixes the near horizon behaviour of the 
perturbation. But this is all that is necessary to calculate the canonical momentum, as we've 
seen in the previous sections. For a two derivative theory 

1 V 

TTo; = V^g''''dzcl)u,{z) = iiujcpo) . (6.6) 

K L07rG]Sf 

The existence of a double pole hasn't changed anything. The canonical momentum is still 
regular at the horizon, and from it we can read off the transport coefficient associated with 
(j), ^ = V/k. This immediately implies the universality of r//s = l/47r in these theories. 

It is clear that this type reasoning does not change when we go to higher derivative 
theories, and the generalized canonical momentum method described in section 3 is still valid 
for these backgrounds. 

6.1 Shear viscosity 

It seems that in these extremal backgrounds the analytic formula for the shear viscosity and 
the pole method have to be modified. This is because their deduction crucially relied on the 
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presence of a simple pole in the action, which is no longer the case here, since the Il{z)(f)'{z) 
term is now 0{l/z'^). However the structure of the canonical momentum as a function of (j)'{z) 
is necessarily the same at zero and finite temperature. That is, at the horizon the canonical 
momentum is always of the form 



Uiz) = - v^5-5.<A(^) = ^V^5^%/-^<^o. (6.7) 
K K \l gtt 

Whether Qzz has a simple or a double pole is irrelevant in the above calculation, since in the 
end we are only interested in finding the coefficient k. This means that the passage from 
finite to zero temperature is smooth, and so we are free to compute the shear viscosity at 
finite temperature and then take the zero temperature limit. This will necessarily yield the 
correct strictly zero temperature result. Therefore we can simply take the background off 
extremality, do the computation there and then extrapolate to T = 0. 

While this suffices to make the analytic formula and the pole method work, we can 
actually do better, since we are not fully exploring the symmetries of the problem in this 
fashion. Since we are only interested in the near-horizon behaviour, we can directly use the 
radial formula for shear visocisty in the AdS2 x R'^~^ geometry. It is simple to see that this 
leads to 

V / a^* \ 

V = -n-FT- ^X^lv - X^l, - 8— . (6.8) 



The coefficient (9^a^* has vanished as there can be no z dependence by symmetry. 

It should be emphasized that this formula is only valid assuming that there are no covari- 
ant derivatives of the gauge field strength in the action; otherwise the shear viscosity receives 
extra contributions from this type of terms. As for the pole method, it can also be directly 
applied in the near horizon geometry, and it works for any higher derivative lagrangian. We 
do need to work at finite temperature in order to have simple poles. This is accomplished 
with the AdS2 black hole geometry, 

ds^ = -vi -{z^ - zl)dT'^ + — 2 + V2{dxidx'-), (6.9) 

V ^ ~ ^0 / 

Fzr = Q, (6.10) 



which is obtained by taking the T ^ and near horizon limit simultaneously as in [|29| . 

Higher derivative corrections modify vi,V2 but leave the form of the metric invariant. 
Then the residue formulas are directly appliable, and the final result does not depend on the 
parameter zq, as expected. 

Let us see how this works in practice. Consider the action 

1 /■ , 12 1 



IGttGtv 



I d'xV^ + - ^2 + ^L^RabcdR"'"'^ (6.11) 



with 7 <C 1. This action has a solution of the type AdS2 x R^, corresponding to the near 
horizon limit of the extremal AdS-Reissner-Nordstrom black hole. The parameters vi,V2,Q 
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are given by 



^1 = ^, ^^2 = ^, g = ^(1 -247) (6.12) 



The temperature is then given by T = zq/{2'k). Adding a shear mode perturbation of the 
form </>(z) = e~*'^* leads to a pole of the lagrangian at z = zq. The residue is 

Re=..„41_.„ = -I^J^ Cif (1 - 487). (6.13) 
Applying the time formula ([4.5D gives 



L_m'(l_487). (6.14) 



This result agrees with the one obtained in [25|. Alternatively, one may use the formula ( [6.81 ). 
The various quantities are determined to be 

2X^%-X^ly = l, «^* = ^, (6.15) 

leading to the same result. We have tested both these methods for various corrections by 
comparing them with the generalized canonical momentum method, finding always an agree- 
ment. 



6.2 Conductivity 

Now let us consider the conductivity. Here the story is more interesting due to the presence 
of a background gauge field. Because of this the computation of the conductivity is no longer 
determined solely by the horizon behaviour. The background charge acts as an effective mass 
for the gauge field perturbation, leading to a non-trivial flow of the associated canonical 
momentum from the horizon to the boundary. 

In the extremal T = limit, we will still able to unambiguously determine the low 
frequency limit of the conductivity. This is because the frequency dependence of the imaginary 
part of the Green's function is completely fixed by the near horizon behaviour. The fiow 
modifies the norm of the perturbation from the horizon to the boundary, but it cannot affect 
the frequency dependence. As shown in [ p9[] , 

ImGRiuj) = Gr{lo = {))dolmG{Lo) oc 0;^^'= (6.16) 

where do is some frequency independent fiow factor, Q{ijo) is the two point function computed 
in the IR CFT, and is related to the conformal dimension of the operator under question. 
Let us see how this works in practice. 

Consider the Einstein-Maxwell action with a cosmological constant va. d+ 1-dimensions, 
and focus on the AdS2 x R'^~^ solution supported by fiux described in the previous section. 
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Now turn on small perturbations corresponding to the shear channel: 

a,(t,z) = A{z)e-'^^ 

h=[{t,z) = H{z)e^'^^ 

K{t,z) = R{z)e~'^^ 

h^it, z) = (t){z)e-'''^ (6.17) 

The (j){z) perturbation decouples from the others. Working in radial gauge where R{z) = 
we get two equations of motion plus a gauge constraint 

QA'{z) + V2H"iz) = 0, 
QA{z)+V2H'{z) = 0. (6.18) 

Solving the gauge constraint and plugging it into the first equation of motion we obtain 

A"{z) + -^A'{z) + ^^-Zl^^^A{z). (6.19) 

This is the equation of motion for a massive scalar field with = Q^/(fiL^). In the 
Einstein-Maxwell case this gives us = corresponding to a conformal dimension 

5 = 1/2 + 1^ = 2. This equation is easily solved, with solution 

A{z) = e^p{iLo/z){z - iuj) (6.20) 

where we have chosen infalling boundary conditions. Expanding around the boundary this 
gives 

A{z)^z+'^ + \^ + ... (6.21) 
2z S 



This implies [^], |3C] that the IR Green's function is given by 

Gu. = iiv^. (6.22) 

The flow equation ( |6.16 ) then fixes 

ImG/?,(w) oc ioo^ Re(cr) oc uj^. (6.23) 

This result is universal as long as the geometry contains a flux supported AdS2 factor. Since 
we are focusing on the imaginary part of Gji our calculations say nothing about the imaginary 
part of the conductivity. Indeed, generally there one expects a pole at a; = 0, which by 
Kramers-Kronig relations indicates a delta function for the DC conductivity [El]. ^ 



®We thank Sean HartnoU for pointing this out. Throughout this section we focus on the low frequency hmit 
of the real part of the conductivity. 
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We have considered taking the zero temperature Hmit with oo/T — > cxd, cj/^ — > 0. The 
analogous computation with uj/T 0, tu/fj, corresponds to taking the zero temperature 
hmit of the finite temperature conductivity. In that case we also get a zero conductivity limit, 
except for the special case d = 2 + 1, where it is a constant [p^ ]. 

However, there is yet another way of taking the zero temperature limit, and that is 
by requiring lv/T ^ zq, uj/ fi ^ 0. In this limit we are working with the AdS2 black hole 
background. ^ 

The equation of motion for the gauge perturbation is then 

^T-) + ^A-(.) + "'-'-''7°"3>- ^(.). (6.24) 

Z Zq z 
The solution of this equation that is infalhng at the horizon is 

(6.25) 



AW = .„exp(--jr/2--|£p 



lUJ z 



with £p an associated Legendre polynomial. This solution has the boundary expansion 

3 iu^ 
2^ Sl' 

from which we read off 



A{z) ~ z + ^ + (6.26) 



Gr{uj) = liu:\ (6.27) 

o 

Therefore in this case we still obtain zero conductivity in the low frequency limit. 
6.3 Conductivity with arbitrary higher derivative corrections 

Now let us see what happens when one includes higher derivative corrections to the Einstein- 
Maxwell action. It is instructive to consider the simple case where the action is defined by 
( |6.11 ). We take the background geometry as in ( |6.1Cl|) ,( |6rT2 ). The equations of motion for the 
perturbations ( |6.17|) can be derived using an effective action approach (see e.g. 0]). Here 
we show the gauge perturbation equation together with the gauge constraint: 

A"{z) + -A'^{z) + — A{z) + 2V^V2^^ = iSa 

z z 

1=A{z)+V2H\z)=^Sb. 

The sources are 

Z^ 

^ ^ -V^z^Ajz) + 12V2 [{u^' - 3z')H'iz) + z'i2H"iz) + zH"'iz))] ^ 



^The other two limits then correspond to taking zo — > where we get back pure ^4^52, or zo —* +oo where 
the AdS2 space disappears and we are left with a regular non-extremal black hole. 
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Working perturbatively in 7, we apply the lowest order equations of motion on the sources 
to simplify them. In the end we get two simpler equations: 



Miraculously, the 7 dependence has completely dropped out. One could in principle expect 
that the "mass term" could receive 7 corrections. The fact that it doesn't means that the 
conformal dimension associated to A{z) does not get renormalized by the curvature squared 
correction to the lagrangian. 

This result seems to be completely general. We have repeated this procedure for a variety 
of higher derivative theories, including terms of type RFF, (VF)^, (VVF)^, i?(VF)^, F^,etc. 
We have also done it in the AdS2 black hole geometry, and have even performed these pertur- 
bative computations to second order*^. The equation of motion for A(z) stubbornly remains 
the same, although the gauge constraint receives corrections at all orders. 

Actually, this result also applies to the shear viscosity calculation. Higher derivative cor- 
rections to the equation of motion for </)(z) defined in ( |6.17| ) can also be treated perturbatively, 
and there one also gets back the original lowest order equation of motion. Of course there we 
have an understanding on why this is the case, since the equation of motion for a massless 
scalar is essentially fixed by demanding horizon regularity, or equivalently it is fixed by the 
conformal symmetry of the background. Even though the same argument necessarily fixes 
the form of the equation of motion for the gauge perturbation A{z) deep inside the AdS2 
bulk, the situation is different because of the presence of a background charge. This acts as 
an effective mass, producing a flow from the AdS2 horizon to the boundary. There seems to 
be in principle no reason why higher derivative corrections shouldn't modify this flow. 

Even though we haven't been able to formally prove that the equation ( |6.29| ) is unaffected 
by perturbative corrections to the lagrangian, our results strongly indicate this is the case. 
We therefore conjecture that this is true to all orders in perturbation theory. This means that 
from the point of view of the IR CFT, the scalar operator ax{t,r) has protected conformal 
dimension 6 = 2. 

This result has important consequences for the conductivity of field theories at extremal- 
ity. Since the frequency dependence of the Green's function doesn't change, we conclude that 
the real part of the conductivity scales like in the low frequency limit, to all orders in the 
higher derivative expansion. 

^Notebooks are available from the author upon request. 




1 + 247 



A{z)+V2H'iz) = 0. 



Plugging the second equation into the first we get 




(6.29) 
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7. Discussion 



In this paper we have shown that the approach of ||21[ is generahzable to higher derivative the- 
ories. Thermal retarded Green's functions are given by the generahzed canonical momentum. 
In the hydrodynamic, low frequency limit where the flow of this quantity is trivial, regularity 
at the horizon completely determines its form up to a constant. Similarly, we've also seen 
how at the horizon the action is also specified up to that same constant factor, the effective 
horizon or membrane coupling, which essentially determines the associated transport coeffi- 
cient. We would like to note however, that our method only applies for first order transport 
coefficients, that is, those that can be obtained by a Kubo formula in the low frequency limit. 
Higher order coefficients require going beyond the low frequency and zero momentum limit, 
and in this case one has a non-trivial flow from the horizon to the boundary which must be 
taken into account. 

To determine the membrane coupling one exploits the form of the action at the horizon to 
generate poles in the lagrangian density. The residue of the simple pole then simply determines 
the associated transport coefficient. The pole method given by equations (4.4) and ( [4.5| ) is 
generally applicable to any higher derivative theory, and requires only the evaluation of the 
lagrangian on a perturbed background, making it extremely efficient. We hope this method 
will facilitate future computations of corrections to transport coefficients. 

We have also derived analytic formulae for the shear viscosity and conductivity, in equa- 
tions ( p.7D ,(|5l8D,( 5.Il| ). Though these formulae are not fully covariant, it should be possible 
to make them so with some work. It is easily checked that they do respect radial coordinate 
reparameterizations. Our formula for the shear viscosity shows that the proposal of is 
incomplete, as it does not properly account for the contribution of higher derivative terms in 
the action. 

Ultimately our methods work because of the specific constraints that arise at the horizon. 
This seems to be related to an emerging conformal symmetry which completely fixes the 
behaviour of any perturbations. It would be interesting to understand this in more detail. 

We have shown that the pole method for computing the shear viscosity is still applicable 
in extremal black hole backgrounds. Alternatively one can use the simple formula ( |6.8| ). We 
have found that in the dual IR CFT there is an operator which seems to have protected 
conformal dimension to all orders in perturbation theory. It might be possible to obtain a 
proof of this statement, and understand generically which operators satisfy this property by 
performing a dimensional reduction onto an effective two dimensional theory. This would be 
important as in our case this fact is responsible for fixing the low frequency scaling of Re(fT) 
to be 0{uj'^). 

It seems likely that our methods can be extended more or less straightforwardly to back- 
grounds dual to field theories with non-relativistic symmetries j3^,^]. This is because these 
backgrounds have similar horizon structure as in the relativistic case. Also it would be inter- 
esting to see if we can adapt our methods to study fermionic correlation functions. We leave 
these and other questions for future work. 
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A. Derivation of the formulae for the shear viscosity. 



We start with the geometry ( |2.lD . The Christoffel symbols are 

= + = 1 + f'{z) = p'{z). (A.i) 

The strategy is to find the simple pole contained in the effective action for the shear mode 
( |5.5D . These come from ^'^ factors upon plugging in the near horizon form of the perturbation: 

V^gg^^cPU^)cP'_Jz) - (A.2) 

It is also useful to notice that at the horizon we have 

VzF^ ~ -VtF\0{l/z), V,,F-; ~ 0{1) (A.3) 

Start by considering the two derivative terms. We define X2 = y,xi = x. With Ay{z) = 
(j){z) we obtain 

^2pj^abcdp^^p^^ = 4X^%g'^cl,'Jz)^'_Jz) 
e'PX^'^nan,e'PF,,F/ = (X^ , + X^y,y)g^^cP'J^)ct>'_Jz). (A.4) 

Therefore the two derivative part of the effective action leads to 

^{-2X^y^y + X^ly) (A.5) 

where we have used X^^^^ = X^^yz by rotational invariance. Next up are the four derivative 
terms. These are 

+ 2Y^^^ye-^PVmFnyV-Ae^'Fzy) + Y^^:^yV ^FnyV^F'^y] . (A.6) 

Here m,n,p,q indices refer to all coordinates with the exception of z. For brevity we will 
now define Y"'^ = Y^il^^y d'^Qaa- Taking the near horizon limit and keeping only the poles we 
obtain: 

^ = 6/9'(y^^ - Y^'){v zF^y)F\ + (y** - Y^'){VtF\f 
+ 2(y*- - y-)(VtF* )(v,n) - 2(y^-» - y*-0(VtFM(v,,no. (A.7) 
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The differences of Y tensors correspond to the a parameters defined in (5.1). In the above 
expression the second and third terms are actually double poles. This imposes the regularity 
condition y** = Y^^. The other term is harmless as the double pole arises from the presence 
of a total derivative. Here we simply take the simple pole part, as this is all our method 
requires. It is easy to show that 

V,F; = -^e-2f (1 + 2zg' + 2zp') </)'(z) 
VtF\ = ^e'^3 {l + 2zf) <p'iz) 

The single pole part is extracted as follows: 

^a'\VtF%)iV,Fl) ^ Oil/z') + \^ge-^3 [d,a^' - a^'e^m] if^^'{z)<p'iz)) (A.8) 



where we have used 



lim R = -e-2s(3/' - g' + 2{d - l)p') 



. Putting all the ingredients together the effective two derivative lagrangian reads 



r(i) — ^ ^ ( nv^y — x^y — Ap-'^9o 
^ ~ ustt^tgn z [ ^ 



' iA.9) 



Using ( p~4| ) then leads to the result ( |5.7| ), except for the a* terms, defined as 

a* = y^^'-y*^\ (A.IO) 



These are actually zero by regularity. This follows from deducing the time formula ( p. 81) , which 
can be done in an analogous way to what we have just done. However, in this deduction the 
a* coefficients do not enter. Equality of the two formulae then implies the constraints: 

yxz ,yz _ yxz ,yt _ yxt ,yz , yxt ,yt 

zy zx zy tx zy tx ' zy zxi 

yxz ,yz _ yxt ,yt _ p 
zy zx ^ ty tx ^ 

yxxi ,yz _ yXXi ,yt _ p 
Xiy zx ^ Xiy tx ^ 

X\-X^% = 0. (A.ll) 

The equalities above should be true for any background and lagrangian, sufficiently close to 
the horizon. 
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B. Functions in 4.12 



The functions in equation ( [4.12 ) are given by 

A{z) 



B{z) 
C{z) 



-gg 



D(z) 



-g{gzz)2 
zz\3 



-g{g'') 



1 - 647(18 - 58z + 77^2 - 48z^ + 12z^) 
647(1 - z)(24 - 1062 + 1772^ - 124z^ + SW^) 
-1927(1 - zf{A - 22z + 472^ - 36z^ + 9z^) 
967(1 - zf 

-327(1 -z)(3- 10z + 5z^) 
-47 
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